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Dynamics of the free jets from nozzles of complex geometries 
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Abstract 

The dynamics of the coherent structures in jets generated by nozzles of different shapes is analyzed through 
DNS at Rcd^ = 565, by considering circular, square, fractal and star-like nozzles. The jets generated from 
orifices with corners, undergone a rotation proportional to the corner angular width: 6rotation = Ocomer/'^- The 
velocity at which this rotation occurs is also affected by the angle of the corners, being faster for fractal 
and star-like nozzles which have small 6comer ^ Therefore it has been found that the velocity of the rotation 
is associated with enhanced spreading and entraining characteristics. The jet evolution and its rotation are 
dictated by the vorticity field and, in particular, by the positive and negative cOx layers generated at each 
corner. The comparison between the fractal and the star-like jets at this Re, suggested that the effect of the 
smaller scales generated by the fractal nozzle does not play a role in the development of the jet, that evolves 
as the star-like one. 


1 Introduction 

The starting jet formed by the sudden discharge of hot reaction products into a fuel-air mixture may have 
a large influence on the combustion process [1], with important implications in the transport, handling and 
storage of fuels, particularly hydrogen, and may also find application in engine ignition systems [2]. The 
control of the jet evolution is strongly dependent on understanding the dynamics of the vortical structures, 
because the spreading is affected by the formation, interaction, merging, and breakdown of these structures. 

The starting jet generated by a circular orifice shows a leading vortex ring followed by a slender jet 
stem, and the associated flow dynamics has been studied for constant density configurations by [3]. They 
observed that, as the boundary layer separates at the orifice, the vortex sheet rolls up to form a toroidal 
vortex that travels downstream, entraining the outer fluid. As a result of the roll-up, a mixed core of reactants 
and combustion products should form at the jet head, providing a precursor ignition kernel where chemical 
reactions are enabled by the high temperature. Successively, [4] and [5] found that the entrainment of jets can 
be enhanced by using non-circular nozzles. The reason for the enhanced mixing and entrainment properties 
has been attributed to the characteristic rotation of the jet cross section, that occurs during its development in 
the stream-wise direction. This axis rotation results from self-induced deformation of vortex rings with non- 
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uniform azimuthal curvature, and it has been observed using elliptic nozzles and nozzles with corners (16], 
17]). In particular, laboratory experiments by 16] have shown that such rotation in square jets is responsible 
for enhancing the mixing of fluid in the neighborhood of the corner regions and near the jet outlet. 

In the present work, DNS simulations at low Reynolds number have been performed to study the near- 
held evolution of the vortical structures generated by nozzles of complex shapes. 


2 Numerical experiments 


The details of the numerical scheme to solve the Navier-Stokes equations, together with the immersed bound¬ 
ary method used to reproduce the interaction between the how and the solid, are described in 18]. The shape 
of the orihces considered are given in hgure 1. The fractal orihce has been obtained with three iterations of 
the basic triangular geometry and all the orihces have been designed to have the same area, thus the solidity is 
the same in all the cases. De is the equivalent diameter of the jets and it is De = 2.257 since the area is equal to 
4. At the inlet an uniform velocity prohle is imposed and the orihce is located at = 2.25. The simulations 
were performed at Reo^ = UooDelv = 565 (Uoo = 1) in a computational domain Li xL 2 xL 3 = 37rx27ix27i 
discretized with n\ x 112 xn^ = 3^5x251 x 251 points. The orihce is located in y - z and v is the downstream 
direction. Periodicity is assumed in y and z. In the plots the v-coordinate is normalized with De starting from 
X^\ V* = {x- X]^)IDe. 



Fig. 1: Geometry of the orihces {a) Js, (b) Jc, (c) Jf, (d) Je and 5 '- and /-directions. 


3 Results 

These DNS allow to study the vortex dynamics governing the evolution of the jets. At Re = 565, the how 
is laminar at the exit from the nozzle, and therefore it is possible to understand the causes leading to the jet 
deformation. In hgure 2 the surface contour ofU= Uoo/2 are given for the square, circular, fractal and star- 
like jets in yellow. The how generated by the circular orihce (Jc) is spread by viscous diffusion and maintains 
its circular shape downstream; in the other cases the shape is dictated by the hole only in the proximity of 
the orihce, and is lost is a short distance. At the end of the domain the how is characterized by a number of 
corrugations that are rotated with respect to their initial position. For Js these are rotated of 45°; for Jf only 
the six corrugations, related to the largest scales, are recognizable from the contour of U, and these are rotated 
about 30° with respect to the initial conhguration. The same rotation occurs for J^. Therefore in the cases 
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with corners, the rotation of the jet is proportional to the number of spikes, in particular: Orotation = Ocomerl'^- 
This rotation was studied experimentally and numerically by [9] and [10]; in particular [10] considered 
square, rectangular and equilateral triangular orifices, observing a 45° rotation for the first, an axis-switching 
(90° rotation) for the second, and an overturning of 180° for the latest. Considering only the initial and final 
position, the overturning of the equilateral triangle is indeed equivalent to a 30° rotation, with Orotation = 
0corner I Therefore also the work of [10] corroborates the present results of Orotation = Ocomerl'^, which is 
valid only for the shapes with aspect ratio equal to 1. In fact the rectangular and the elliptic orifice, undergo 
a different Orotation- 


(a) (b) 



Fig. 2: 3D Contour plot of U = f/oo/2 for {a) Js, (b) Jc, (c) Jf and (J) in yellow, overlapped with 
= +0.5 in red and blue respectively. 0 < x* < 3.18. 


The location where the rotation of the jet terminates is also affected by the shape, and this length can 
be determined by monitoring the stream-wise variation of the half velocity width ri /2 measured on the 
smallest (^'-direction) and largest (/-direction) radii defined in figure 1. In figure 3a we reported the variation 
of the half-widths ^ 1/2 non-dimensionalized by De for Jg, Jf and while for Jc we reported 

only because the jet remains circular. The half-width ri /2 of the jet, at the given stream-wise location, 
is the distance from the center-line at which the axial velocity drops to half of its center-line value UcU If 
becomes larger than means that the jet has inverted its long and short diameters, and the rotation 

Orotation = Ocornerl'^ is Completed, becomes larger than at x* =0.513 for Jg, x* = 0.260 for Jf and 
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X* = 0.275 for thus Jp and complete their rotation before Jg . “ ^ 1/2 ^ fundamental quantity 

because indicates the stretching of the corner portions of the jet and, being approximately equal for all 
the jets, where grows faster, the corners undergo a stronger stretching. After the initial region where Jp 
is spread faster than /^, at x* 1.5, roaches its maximum, while it continues to grow for Jg . At the end 
of the domain of Jg is larger than in Jp because the former maintains its coherent shape longer, and the 
spreading is slower; on the other hand Jp, as it will be shown later, produces an intricate flow pattern, where 
the vortex interactions make the jet to lose its coherent shape, and decreases. 

To understand whether small scales corrugations play a role at small Re the fractal jet has been compared 
with / 6 - Figure 3a does not show large differences, suggesting that the smallest scales generated by the 
fractal jet die and do not influence the evolution far from the nozzle. The simulation allows to attribute 
the initial highest spreading of Jp, not to the fractal corrugations, but to the angular width of the largest 
corners. In figure 3a it is reported also ri /2 for Jp to show that with this Reynolds number the radius of the 
jet remains almost unchanged along v*. 



Fig. 3: {a) Evolution of the jet half-width normalized with De vs stream-wise direction v*. Symbols: (■): 

for Jg ; (o); for Jf\ (a): for J6\ (><)• ri /2 for Jp. Dashed lines for the (^) Evolution of {Qq) 

and {Qx) vs the stream-wise direction v*. Symbols: (□): {Qx) for Jg ; (■): {Qq) for Jg ; (o): {Qx) for Jp', (•): 
{Qe) for Jp', (O): {Qx) for /6;(^): {Qe) for Je\ (A): {Qx) for Jc', (a): {Qq) for Jp. 

In the following paragraph the causes leading the jets to have different spreading are explained by looking 
at the vorticity fields which, being affected by the corners, dictate the evolution of the velocity structures pre¬ 
viously described. The effect of corners is enlightened in figure 2, where the positive and negative contours 
of (jJx are respectively coloured in red and blue, and are overlapped to the contour ofU = Uool2. The figure 
clearly shows that in Jg, Jp and J^, the thin layers of ojx produced at the corners, are responsible for the 
rotation and successive stretching of the jets. The formation of oJx can not be instead observed in Jp. At the 
nozzle exit, in correspondence of every corner, a pair of thin layers of ±cl>x vortices forms. While these are 
convected downstream, the positive vortices come up against the negative ones generated by another corner, 
thus forming a new couple. This interaction between ojx of opposite sign stretches the spikes and spreads 
the jet. Eor the fractal Jp, ±(jl)x of different size are produced at the inlet but, due to the low Re, the small 
patches die and do not affect the flow; therefore at a certain distance the jet is affected only by the six largest 
structures and does not differ from J^. 

In addition to these vortical structures, at the end of the computational domain of Jp (figure 2c) and 
/6 (figure 2 d), other structures with ±cl>x, that are not connected to the layers generated at the corners, are 
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visible. These structures are due to the vortex dynamics that, through the processes of vortex-stretching 
and vortex-tilting, generates ojx from the other vorticity components. In fact at the inlet ojy and are also 
generated by the orifice, forming a corrugated sheet of intense vorticity. Such structures for the three jets 

Jc, Js and Jp are shown in figure 4, where isocontours of Qq = and Qx = in the plane yz are 

plotted at two distances v*. The quantity Qq for Jc shows a structure with the same shape of the vorticity ojq 
for a circular jet, evaluated in cylindrical coordinates [11]. In the figure the region with a clustering of the 
contour lines are those with the highest values of and represents the boundaries of the vortex sheet. The 
increase of the width of this region with x* indicates the jet Jc is spreading both towards the center-line and 
the surroundings. The other jets present a more complicated pattern at x* = 3.18, because the jet spreading 
is greater than for Jc, and the fiuid transport between the jet and its surroundings is enhanced. This fast 
mechanism, triggered by the Qx formed at the corners of the nozzle, is the cause for the intricate vorticity 
pattern, where also Qx not generated at the inlet by the orifice is present. 



Fig. 4: Contour plot of Qq = + normalized with its maximum value in the range [0.1 : 1] with 10 

contour lines in red for {a) Js , (b) Jc and (c) at x* = 0; (J), (e) and (/) at x* = 3.18. The plot is overlapped 
with the contours of Qx normalized with its maximum value in the range [0.3 : 1] with 7 contour lines in 
blue. 


Figure 3b reports the evolution of (Qq) and (Qx) averaged in the plane yz and in time, along the stream 
wise direction. The values achieved by (Qq) are larger than (Qx) for any jet, and the plot shows that (Qq) oc 
^I ncomer', in fact the highest values are achieved by Jp and J^ that have a smaller 6comer with respect to 
Js . The plot shows that {Qx) is null along the whole domain for the circular jet. In the other cases {Qx) 
is approximately constant, until x* 2.25. The subsequent growth is due to the vorticity dynamics that 
produces {Qx) from {Qq) by vortex stretching and vortex tilting. The high {Qq) is characterized by a first 
region where this quantity decreases, followed by a second region where it grows linearly. At the end the 
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decrease is mainly due to the viscosity. The slope and the length of the linear region depend on the geometry 
of the orifice, in fact aj^ = 0.14 and ajj, = 0.14, while aj^ is 0.10 and the value of {Qq) for Jc remains 
approximately constant = 0.01). 


4 Conclusion 

The present DNS of the jets produced by nozzles of complex geometries at Rcd^ =565 allowed to visualize 
the complex vortex dynamics, to increase the understanding of the jet evolution, in order to enhance the 
entrainment and the mixing. The jets generated from nozzles with sharp corners evolve differently from 
circular jets, that maintain an axisymmetric circular shape, due to the axis switching of the jet. This rotation 
is produced at the nozzle and is initiated by the vorticity self-induction process. The rotation of the jet has 
been related to the angular width of the corners {Orotation = Ocomerl'^) and, from the data in literature, the 
relationship between Orotation and Ocomer seems to be valid for nozzles with aspect ratio 1. The velocity at 
which the rotation occurs can be measured by the growth rates of the 5'- and /-diagonals and it has been 
shown that, while similar for all the nozzles with comers, strongly affected by the shape. The 

growth of associated with the stretching of the corrugations, in fact in Jp and grows faster, the 

corners undergo a faster stretching, and this process contributes to enhance the spreading. The deformation 
of the jet is linked to the vortex dynamics because patches of ±0Jx are formed at the corners. These thin layers 
are convected downstream and in turn deform the jet. The spikes of the jet are stretched by these vortices 
and consequently the structure is spread more rapidly, higher is the strength of the ojx patches. As a result 
the fluid transport between the jet and its surroundings is enhanced, so as the mixedness. The low Reynolds 
number assumption guarantees that the flow remains laminar in the near-field and, in such condition, the 
wider range of scales generated by the fractal nozzle, does not affect the flow, therefore the evolution of Jp 
and /6 does not show large differences. 
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